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We study the motion of small particles in a random turbulent flow assuming linear law of friction. 
We derive a symmetry relation obeyed by the large deviations of the finite time Lyapunov exponents 
in the phase space. The relation applies when either the statistics of the strain matrix is invariant 
under the transposition or when it is time-reversible. We show that, as a result, the Lyapunov 
exponents come in pairs which sum is equal to minus the inverse relaxation time of the particles. 
We use the pairing to consider the Kaplan- Yorke dimension of the particles attractor in the phase 
space. In particular, the results apply to case of the fiow which is white noise in time. 



PACS numbers: 05.40.-a, 47.55. Kf, 47.52. +j, 47.27.-i 

Recently much interest was raised by the Gallavotti- 
Cohen type fluctuation relations [l[ for the dissipative 
dynamics, see e. g. [H and references therein. The rela- 
tions describe the large deviations of different observables 
and apply to the chaotic dynamics with time-reversal 
symmetry Q. An example is the relation on the en- 
tropy function describing the large deviations of the finite 
time Lyapunov exponents [3, 4]. It is natural to inquire 
whether other symmetries of the dynamics also generate 
the corresponding fluctuation relations. Here we answer 
this question affirmatively in the context of the chaotic 
dynamics describing the motion of inertial particles in 
a turbulent flow. We show that two different statisti- 
cal symmetries of the strain matrix of the flow lead to 
a new fluctuation relation for the finite time Lyapunov 
exponents. This relation implies that the spectrum of 
the Lyapunov exponents has a reflection point: the ex- 
ponents come in pairs which sum is equal to minus the 
inverse relaxation time of the particles. 

The behavior of small inertial particles in turbulence 
is now enjoying much attention from the researchers, see 
e. g. 0, a Hi, S M El, [H, [13, [II- The problem 
has many applications including rain formation in clouds 
0, El , ocean physics and engineering Q . Theoretical 
study of the problem mostly involved modeling turbu- 
lence by a white noise in time, Gaussian velocity fleld, 
known as Kraichnan model [j, 0]- Even in that case 
theoretical study is rather difficult and analytic results 
were mainly obtained for the one-dimensional case 
In particular, the spectrum of the Lyapunov exponents 
is unknown in dimension larger than one and probably 
cannot be found analytically 11 ill . In such a situation nu- 
merical studies are necessary [ill [l^ . The pairing of the 
Lyapunov exponents proved here for a class of flows that 
includes Kraichnan model is expected to be a useful tool 
in such studies. 

Let us briefly review the study of the pairing of the 
Lyapunov exponents for dynamical systems. The pairing 
is a well-known property of the Hamiltonian systems: the 
exponents come in pairs with zero sum, see e. g. [l5| . 



This is a consequence of the symplectic property of the 
linearization of the flow, see [15| and below. The first 
observation of the pairing for the dissipative dynamics is 
due to U. Dressier ^16^. It was shown in jT6l| that the 
addition of the linear friction term to the Hamiltonian 
dynamics merely shifts the (zero) reflection point of the 
spectrum by minus half the friction coefficient. Later, 
the relation appeared in the context of the thermostatted 
systems which are Hamiltonian systems with an added 
non-linear friction term 17, 18|- When the friction term 



is such as to make the kinetic energy of the system con- 
stant (the so-called isokinetic thermostat |18l| ), i t was 
shown by C. P. Dettmann and G. P. Morriss in [19| that, 
after the exclusion of two zero Lyapunov exponents, the 
remaining Lyapunov exponents obey the conjugate pair- 
ing rule. The exponents come in pairs which sum to 
the finite-time average of the phase space contraction. 
Both Dressier and Dettmann-Morriss relations apply for 
the finite time Lyapunov exponents and express the fact 
that Jacobi matrix of the flow is conformally symplectic - 
it transforms simply under the similarity transformation 
by the symplectic tensor, see [l^, [l^ and below. On the 
other hand, it was argued that for the isoenergetic ther- 
mostat that conserves the total rather than the kinetic 
energy, the exact pairing of the flnite time Lyapunov ex- 
ponents should not hold [2^. Still, the pairing may be 
true in the limit of large times (it is almost certainly true 
for both large times and large number of particles when 
equivalence between the isokinetic and isoenergetic ther- 
mostats is expected) (i^l . The isoenergetic thermostat is 
an example of the situation where Jacobi matrix is not 
conformally symplectic and yet the pairing is expected to 
hold asymptotically. To treat such cases different tech- 
niques are needed. For the case of a thermostatted sys- 
tem where the interactions can be described by white 
noise, the use of Fokker-Planck equation allowed to prove 
the pairing rule in [2l|. Here we suggest a different ap- 
proach to pairing in different situations that does not 
involve the assumption of the white noise. We note that 
to prove the conjugate pairing rule for the infinite time 
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Lyapunov exponents it is sufficient to show that Jacobi 
matrix is only statistically conformally symplectic. Here 
statistics is defined by sampling the initial points with 
the natural measure. Indeed, after it is proved that Ja- 
cobi matrix is statistically conformally symplectic, one 
can use the multiplicative ergodic theorem by Oseledec 
(2^ that ensures that the part of the matrix defining the 
Lyapunov exponents does not fluctuate in the limit of 
the large times. Thus in the limit of the large times the 
statistical property brings a deterministic pairing rule for 
the infinite time Lyapunov exponents. 

We pass to the analysis. We consider the motion of a 
small spherical particle immersed in a fluid. We assume 
the d— dimensional flow of the fluid u{r,t) to be smooth 
below some scale rj (if m is a developed turbulent flow 
then 77 is the Kolmogorov scale). Due to the finite inertia 
the particle deviates from the carrying fiow. Provided 
the particle is sufficiently small and heavy, the drag by 
the flow is well described by the linear. Stokes' law of 
friction, see below and for the validity conditions. 
Designating the particle position and velocity by r{t) and 
v{t), Newton's law reads 



V — u{r{t), t) 



(1) 



Here r = {2 / %){po / p){a? / v) , where po and a are the par- 
ticle density and radius, while p and v are the fluid's den- 
sity and kinematic viscosity. Above the particle radius 
should be much smaller than 77, the Reynolds number as- 
sociated with the particles' motion relative to the flow 
should be small and the particle should be sufficiently 
heavy for effects like added mass to be negligible, see 
[l3i] for the details. The trajectories X{t, x) in the phase 
space of Eq. ([T]) obey 



dX{t,x) 
dt 



V{X{t,x),t), X{Q,x)^x, (2) 



where V = {v, —[v — u{r,t)]/T). Note that the flow in 
the phase space is dissipative: V -V = —d/r. Below we 
shall assume that u{r, t) is a random flow with stationary 
statistics. As a concrete example we shall use Kraichnan 
model where turbulence is modeled by a Gaussian veloc- 
ity field u{r, t) with zero average and pair correlation 

{ui{ri,ti)uj{ri,t2)) = S{t2-ti) [Vo5.,j-K.,j{r2 - ri)] . 

Above the tensor Kij{r) is assumed to have finite second 
derivatives at r = in conformance with the assump- 
tion that the flow is smooth. Due to the randomness of 
u{r, t) the flow in the phase space is chaotic for typical 
realizations and it can be characterized by the spectrum 
of the Lyapunov exponents which we introduce now. We 
consider the Jacobi matrix Wij{t,x) = djXi{t,x) deter- 
mining the evolution of infinitesimal volumes in the phase 
space. Consider as an example the distance {R, V) be- 
tween two infinitesimally close trajectories in the phase 



space that satisfies 

dR 

'dt 



= V, 



dV 
~dt 



V (jR 



(3) 



where Uij — V jUi{r{t),t) is the strain matrix evaluated 
at the trajectory of the particle. The solution can be 
written as [R{t),V{t)] = W{t)[R{0),V{0)]. Jacobi ma- 
trix obeys the corresponding equation 



dW 

~dr 



1 



W, W{0) = I. 



(4) 



where each entry in the stability matrix above is a d x d 
matrix and lij = Sij is the 2ci— dimensional unit ma- 
trix. According to Oseledec theorem the evolution of 
distances has well-defined exponential behavior at large 
times corresponding to the existence of the finite ma- 
trix limit \imt--foo^i^W^{t,x)W{t,x)/2t. Assuming that 
the system is ergodic, the eigenvalues Xi of the limiting 
matrix are the same for almost every realization of the 
velocity and they can be considered independent of x. 
These eigenvalues define the Lyapunov spectrum of the 
system. Conventionally we arrange the exponents in the 
non-decreasing order, Ai > A^+i. We introduce the fi- 
nite time Lyapunov exponents /Lti(t, x) as the eigenvalues 
of [\nW'^{t,x)W{t,x)]/2t, so that Muit-,^ p,{t,x) = 
(again pi > Pi+i). The following constraint follows from 
V -V = -d/r: 



(5) 



Below we study the statistics of W{t) as it is obtained 
by averaging over the statistics of the random strain a in 
Eq. ^ . At large times the statistics of the latter becomes 
stationary and, by ergodicity, independent of x (at small 
times the statistics is generally time-dependent because 
fixing the initial condition at t = breaks the transla- 
tional invariance in time). Since we are interested only 
in the asymptotic statistics of pi{t) at large times, then, 
with no loss, we will use for the statistics of a in Eq. ^ its 
large time, stationary, asymptotic form. The stationary 
statistics holds automatically if we consider the fluctua- 
tions in the steady state when initial points are sampled 
with the natural measure. The simplest example of the 
statistics of a is provided by Kraichnan model where a 
is Gaussian with zero average and the pair correlation 

(CTy(tl)(T,„„(<2) - S{t2 - ii)V,V„i^,™(0). (6) 

Note that when a in Eq. ([5]) is a white noise, still there 
is no need to specify the interpretation convention (say 
Ito or Stratonovich) because the correction term, which 
arises when changing the interpretation convention, van- 
ishes identically, see [lH. We now pass to the deriva- 
tion of the pairing rule. We follow the same lines as in 
0, [3 • introduce the symplectic tensor J by 



J = 



1 

-1 



(7) 
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that obeys J^J = 1. It is easy to see that Ws{t) 
-t/T]j'^W^^''^ (t)J obeys the equation 



exp[ 



1 



dt V2 i 



Ws, WAO)=I- 



(8) 



Above the subscript s stands for symplectic, see be- 
low. Eq. ^ is the same as Eq. ^ up to the 
change a — » <t^ . Noting that the eigenvalues of 
[W^{t)W,{t)]-^ = eM2t/T]JW'^{t)Wit)j'^ coincide 
with those of exp[2t / t]W'^ {t)W (t) , we conclude that the 
spectrum of [In WT (t) Ws (t)] / 2t is given by -fJ-i{t) - l/r, 

cf. [Uliillii. 



The simplest situation is the one of the potential flow, 
u{r, t) — V(l){r, t). In this case Eqs. ([T]) are Hamiltonian 
up to the linear friction term —v/t, the situation con- 
sidered in 16]. We have — a which leads to Ws{t) = 
W{t) and [lTiW^{t)W{t)]/2t = [InWj' {t)Ws{t)]/2t. The 
fact that W differs from J^W~^'^{t)J only by a multi- 
plicative constant is sometimes described by saying that 
W is conformally symplectic It follows from the 

above that the spectrum of [InW'^ {t)W{t)]/2t has a re- 
flection point at — 1/2t and the so-called strong conjugate 
pairing rule [20[ holds: fj,i{t) + fJ-2d+i-i{t) = — I/t. Such 
pairing of the spectrum holds for every realization and 
any time t. Taking the limit of large t we find the pairing 
of the Lyapunov exponents + X2d+i-i = ^^/t- The 
relations agree with Eq. (O . 

For turbulence, where vorticity is present, a'^ ^ <j and 
we cannot apply the above relations. Yet, as we pass 
to show, it is enough for the pairing rule for Xi that W 
is conformally symplectic only statistically. If statistics 
of (T is invariant under transposition cr cr^ then it 
follows from Eqs. (g]) and ^ that W{t) and Ws{t) = 
exp[—t/T]j'^W^^'^{t)J are equal in law. This gives that 



lnW^{t)W{t) \nJ^W~\t)W-^-^{t)J 



2t 



2t 



(9) 



where we designated the equality in law of random vari- 
ables by ~. In the limit t —>■ oo the matrices above have 
finite limit that produces the pairing rule 



Xi + X2d+l-i 



(10) 



At finite t, Eq. ([9]) implies a symmetry of the probability 
density function (PDF) P({^J,i) of n^{t) 



P{fJ-i, ■•,M2d) = P 



1 1, , ^ 
-M2d • (11) 

T T 



We remind that by our convention on ordering of fii the 
domain of definition of P{xi, ..,X2d,t) is xi > .. > X2d- 
Let us rewrite Eq. pip in a form more usual for fluctu- 
ation relations. At large times the probability of a finite 
deviation of jii from their long-time limits Xi decays. It 
was shown in [23| that the decay is exponential when the 



Lyapunov spectrum is non-degenerate, Xi > Xi+i, see 
also [4i| . Accounting for the constraint ^ one has 



d/T] CXp[-tiJ(^i, ..,//2d-l)], (12) 



Here H{xi, X2, ■■, X2d-i) is a convex entropy function 
which is positive everywhere except for Xi = Xi, where 
it vanishes. The quadratic expansion of H near its min- 
imum at Xi — Xi would reproduce the central limit theo- 
rem Q. Using Eq. mi]) wc find the fluctuation relation 

H {ni,..,H2d-i) = H ^^-Ai2d, ■•, ^-/^2^ , (13) 

where /i2d — —d/r — X]i=i ^ l^i- Evaluating the above at 
/ii = Xi one recovers Eq. (fTO|) . 

To discuss the relevance of the symmetry a — > cr^ for 
the turbulent flows, let us first note that since small scale 
turbulence is isotropic, it is natural to consider the ro- 
tationally invariant statistics of a. Then, in Gaussian 
approximation the statistics of a always has the symme- 
try a — » cr^. Indeed the mean value and the pair cor- 
relation function (fully determining the statistics) have 
the general form (cr,:j) = DSij and {(Tij{ti)(Jmn{t2)) = 

A{t2 — ti)6ijS„in + B(t2 — ti)SimSjn + C(<2 — tl)SinSmj 

which is invariant under the substitution a . It 

follows that, in particular, Eq. pUj) holds for Kraichnan 
model with isotropic statistics. 

The fluctuation relation (fTTj) is different from the 
Gallavotti-Cohen relation which for PDF reads as 

P{^il,■■,^i.2d) = P(-Ai2d, ■■, -Mi)exp[X)Mi]: see 0,0]. 
Gallavotti-Cohen relation applies under the condition of 
time-reversibility, which is not obeyed by Eq. ([T]) . Never- 
theless, under the condition that the "generalized time- 
reversal" symmetry a(t) —^ (7{—t) holds, one can derive 
(fTS]) from Gallavotti-Cohen relation in some transformed 
variables. Here "generalized time- reversal" stresses that 
the dynamics in the phase space is not time-reversible 
and we talk only about the symmetry of the strain ma- 
trix. We first rewrite Eq. Q as 



d^R 

dt^ 



1 dR 

T~dt 



aR 

T 



(14) 



We introduce ^p{t) = exp [t/2r] R which satisfies 



df 



1 

4^2 



where the evolution matrix W obeys 





= mW , m 

dt 



1 

473 



It is easy to see from the definition of ip that at large 
times the finite time Lyapunov exponents ili(t) of W{t) 
(the eigenvalues of [InW"^ {t)W{t)]/2t) are simply re- 
lated to ^J,i{t): jli{t) ~ fii{t) + l/2r. Then, at large 
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times, Eq. ((TT|) is equivalent to the ordinary Gallavotti- 
Cohen relation for the PDF of p,i, that is .., fj,2d) = 

P{—fi2d, ^Ai); where we used det W = 1 and ^ jli — 
0. Thus for the entropy function describing the behavior 
at large times, Gallavotti-Cohen relation for W implies 
the validity of Eg. fn]) . 

Let us prove Gallavotti-Cohen relation for W under the 
assumption that the statistics of a{t) is invariant under 
a{t) <y{—t)- We introduce the second time argument 
into W by 



Af + 0(r) and A2d+i- 



l/r - Af + 0(r) for 



dW{t,s) 
di 



rh{t)W{t,s), W,j{s,s)^5,^. (15) 



Writing W{t, s) as the time-ordered exponent 

W = l+ / m{ti)dti+ / dti / dt2rh{ti)m{t2) + .., 

J s J s J s 

it is easy to see that dsW{t,s) — ~W{t, s)m{s). It fol- 
lows that dsW~^{t,s) = m{s)W~^{t, s). Consequently, 
W-'^(t) = W-^(t,0) can be found as M(t = t) where 
drM{T) = ~m{t - t)M(t) and M{t = 0) = /, cf. [1]. 
Introducing the matrix 



N : 



1 

-1/ ' 



(16) 



one finds drN^M{T)N = m{t-T)N'^ M{t)N . The sym- 
metry ij{t) (j{—t) together with the stationarity imply 
the symmetry (t(t) a{t — r) and m{T) —>■ •m{t — r). 
We conclude that the invariance of the statistics of a 
under the time-reversal a{t) —>■ (j{—t) implies the equal- 
ity in law of W{t) and N'^W-'^{t)N. This proves the 
Gallavotti-Cohen relation for W , since the eigenvalues of 
N'^W~^''^ {t)W^'^ {t)N are easily seen to be the same as 
of [WT{t)W{t)]~\ 

Let us repeat that the relation holding as a conse- 
quence of a{t) — > cr[—t) symmetry is Eq. (jl3p and not 
the Gallavotti-Cohen type symmetry for /i^: by itself the 
flow in the phase space is not time reversible. In fact, 
it is hard to formulate the conditions on the statistics of 
tt(r,i) that would ensure the symmetry ij{t) ij{—t). 
One case where the symmetry holds is Kraichnan model. 
Thus the fluctuation relation (fT3)) holds in Kraichnan 
model always, even if the velocity statistics is not ro- 
tationally invariant. 

From now on we consider the implications and the ap- 
plications of the pairing relation for A^ . Let us first study 
the spectrum of the Lyapunov exponents, regarding the 
latter as a function of particles inertia at fixed statistics 
of u{r,t). The inertia can be measured by the dimen- 
sionless Stokes number St = A{'t, where we introduced 

r 7 

the spectrum A^ , i = l,2,..,d of the Lyapunov expo- 
nents of the flow u{r,t). At small St the effect of in- 
ertia is insignificant and the particles simply follow the 
flow. As a consequence of the pairing rule (jlOp we find 



A. . , , ^ 

i — l,2,..,(i. The large negative term — 1/t in A.; with 
i > d describes fast relaxation of the velocity to the local 
velocity of the fiow. As St increases, eventually one enters 
the regime where Air ^ I, which can be called the limit 
of strong inertia, cf. [13] ■ In this limit one can neglect 
the 1/t term in Eq. (|10p and find the simple refiection 
symmetry Xi = —X2d+i-i- Note that the latter approxi- 
mation does not respect the identity J2 •^i — —d/r. 

The above can be used to derive a nice expression for 
Kaplan- Yorke dimension Dky [13 particles at- 

tractor in the phase space. In the limit of large inertia 
one deals with heavy particles which motion is close to 
ballistic. As a result, in this limit the dimension of par- 
ticles' attractor becomes close to the full dimension of 
the phase space 2d. Let us see how the dimension gap 
Cky = 2d — Dky closes. Using the refiection symmetry 



E 



2d-l 



A,, 



Ai > and the identity E = —d/r one 



finds that at large inertia 



Cky « -— , Air>l, 
Air 



(17) 



where we used X2d ~ — Ai. In the case of Kraichnan 
model the scaling argument at large St gives Air cx St^^"^ 
ll|,[il[l3 and we find 



C 



KY 



St 



-1/3 



St 



1/3 



> 1. 



(18) 



We observe that the decay is rather slow. Moreover, the 
above analysis implies that Cky ~ 1 for St^^'^ < 1. As 
a result, one expects the particles attractor in the phase 
space to have a complex structure up to possibly very 
large Stokes numbers satisfying St^^'^ ^ 1, cf. [l^ . 

The above analysis of the limit of the large inertia 
has in fact a universal applicability - Kraichnan model 
becomes valid in the limit of large r for arbitrary flow 
with finite correlations in time |1^ . [l^ . In this limit 
the separation of the time-scale r of the particle and the 
time scales of the flow allows to introduce the effective 
Langevin description in much the same way as it is done 
in the ordinary theory of Brownian motion 25|. Qual- 



itatively, heavy particles perform effective averaging of 
the action of the flow on them, that allows Gaussian ap- 
proximation, cf. [l3l- We make a general conclusion that 
in the limit of large r the codimension of the particles' 
attractor in the phase space scales with r as t~^^^. We 
also conclude that the pairing of the Lyapunov exponents 
holds always in the limit of large r. 

Another range of situations where Kraichnan model, 
and thus the pairing, are of universal applicability is the 
case where a strong external field is present, such as a 
constant gravitational or electric field (in the latter case 
particles must be charged) [3, [3, 27 1. Here the field 
endows the particles with a finite velocity of the drift 
through the fiow (sedimentation) which becomes large at 
large field strength. As a result, the particles cross fast 
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the spatial scale of correlations of the velocity gradient 
and a{t) - velocity gradient measured in particle's frame 
- becomes a fast fluctuating noise describable by Kraich- 
nan model [3, [Mill- 

Finally, let us discuss the applicability of Kraichnan 
model to the case where u{r,t) is a developed turbu- 
lent flow. The limit of large r described above is the 
limit where r is larger than the integral time-scale 
of turbulence and it is not very practical. Yet, it can 
be shown that a modified Kraichnan model holds when 
T is smaller than but larger than the viscous time- 
scale [jJi, i27i . In this version the amplitude of the noise 
appearing in Kraichnan model is a random function of 
time which statistics is not time-reversible. As a result, 
the relations ^TU\i and do not hold in general. How- 
ever, in the limit St^^^ ^ 1 the model simplifies and 
one can use as the amplitude not a random function of 
time but just a random constant [3j[23|- The resulting 
model already has the symmetry a{t) — > a{—t). Since 
for particles in turbulence Air ~ St^''^ at St 3> 1 (but 
T <C ^l) [Hi, then we conclude that = —X2d+i-i holds 
at St^/^ > 1. We also conclude that C'ky ^ St"^/^ at 
St^/^ ^ 1. While the latter limit is rather theoretical 
(St^^® ^ 1 and T <^tL demands very large inertial inter- 
val), it allows to make an important conclusion that for 
particles in turbulence Cky remains of order unit y up 
until very large Stokes numbers, St^^^ 1, see also [ij]. 
Note that here the scaling of Cky with t is not r~^/'^ 
(as claimed above to be universally true in the limit of 
sufficiently large r) because of the limitation t ^t^. 

Let us summarize. We have shown that two statis- 
tical symmetries of the strain matrix, cr and 
a{t) cr(—t) produce a fluctuation relation for the finite 
time Lyapunov exponents which differs from Gallavotti- 
Cohen type relation. In the case of the generalized time- 
reversal symmetry a{t) (t(— i), however, the relation 
follows from Gallavotti-Cohcn relation in transformed 
variables. The fluctuation relation implies that the Lya- 
punov exponents come in pairs which sum to minus in- 
verse relaxation time of the particles (which is formally 
the same result as in [l^). We explore the consequences 
of the pairing of the Lyapunov exponents in different sit- 
uations. In particular, we show a universal result that 
at large inertia the decay of the codimension of the par- 
ticles' attractor in the phase space is slow - the codi- 
mension decays proportionally to t^^^^. For the case of 
developed turbulent flow with very large inertial interval 
an intermediate asymptotic regime may hold where the 
codimension scales as r~^/^. 

We expect that the conjugate pairing rule for the Lya- 
punov exponents of the inertial particles will be useful 
for the numerical measurements of the spectrum of the 
Lyapunov exponents. The question of the existence of 
the fluctuation relations for dynamical systems having 
symmetries other than the time-reversal deserves spe- 



cial attention. We also expect that our approach to the 
derivation of the conjugate pairing rule for the Lyapunov 
exponents will find applications to other situations. 

We are grateful to K. Gawedzki for the kind hospitality 
in Lyon where the collaboration started. We also thank 
him for a remark that helped simplify the exposition of 
the material. 
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